Coherent control of interacting particles using dynamical and Aharonov-Bohm phases 
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A powerful method of manipulating the dynamics of quantum coherent particles is to control the 
phase of their tunneling. We consider a system of two electrons hopping on a quasi one-dimensional 
lattice in the presence of a uniform magnetic field, and study the effect of adding a time-periodic 
driving potential. We show that the dynamical phases produced by the driving can combine with 
the Aharonov-Bohm phases to give precise control of the localization and dynamics of the particles, 
even in the presence of strong particle interactions. 
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Introduction ~ Experimental advances in produc- 
ing low dimensional semiconductor nanostructures has 
given the opportunity of studying quantum transport in 
regimes ranging from single-particle to strongly corre- 
lated. The excellent coherence properties of these de- 
vices, together with the degree of control over their ge- 
ometry and specifications, makes them ideal candidates 
for studying coherent transport, where quantum interfer- 
ence is used to regulate the movement of particles. Such 
control is particularly vital for quantum information ap- 
plications, in which the coherence and entanglement of 
the initial state must be preserved during the evolution 
of the system, and nanostructures represent one of the 
most promising avenues for realizing quantum comput- 
ing elements. 

If we consider a particle hopping on a lattice, inter- 
ference will occur if the hopping acquires a phase factor. 
A direct way of doing this is to apply a magnetic field, 
which produces the well-known Aharonov-Bohm (AB) 
phase. In Refill it was shown that such phases could 
produce a localization effect termed AB caging, in which 
destructive interference bounds the set of sites that can 
be visited by an initially localized wavepacket. For ex- 
ample, in the rhombus-chain lattice shown in Fig[T^, a 
particle initialized on the central site will not spread out 
along the lattice if each plaquette is threaded by a single 
flux quantum ($ = tt). This caging effect has been ob- 
served experimentally in superconducting wire networks 
0, mesoscopic semiconductor lattices [J, and arrays of 
Josephson junctions 4]. The Rashba effect can similarly 
be used to generate tunneling phases, and it has been 
proposed [5| to use the spin-orbit coupling in semicon- 
ductor nanostructures to produce an analogous caging 
effect. 



A contrasting form of localization, also arising from 
quantum interference, is termed ^ "coherent destruction 
of tunneling" (CDT). This arises in systems subjected 




FIG. 1: (a) Schematic arrangement of the "rhombus chain" 
lattice. When each plaquette is threaded by a flux $ = tt, a 
particle initialized on the central site (filled circle) will oscil- 
late to its nearest neighbors (empty circles), but will not prop- 
agate through the lattice. The circled sites form an Aharonov- 
Bohm cage, (b) Occupation probability of the central site 
for a two electron system. The non-interacting case (solid 
black line) displays regular, non-decaying oscillations, indi- 
cating perfect AB-caging. For weak interaction {U — O.IJ, 
dashed red line) the oscillations decay; the caging is no longer 
complete and the particle can leak away. For stronger inter- 
actions {U — J, dot-dashed blue line) the oscillations decay 
even more rapidly, and soon become irregular. 



to a time-periodic driving field. Tunneling processes ac- 
quire phase factors from the interaction of the system 
with the driving, producing an effective renormalization 
of the tunneling. For the most common case of sinu- 
soidal driving, for example, the single-particle tunneling 
is renormalized by the zeroth Bessel function j7o [7j|. 

AB caging is resistant to small quantities of disorder 
Q, but is swiftly destroyed by interactions (see FiglTjD) 
due to the formation of spatially-extended states [q] . In 
this paper we consider a system of two interacting elec- 
trons, and show that by combining a high-frequency driv- 
ing potential with the magnetic flux, stable AB caging 
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can be restored. This occurs when CDT causes the (re- 
pulsive) electrons to bind together into a composite ob- 
ject of charge 2e termed a "doublon" , which can then be 
caged by the magnetic flux. We then go on to consider 
the effect of a low frequency driving field, and show that 
this gives rise to an unusual form of propagation in which 
the doublon moves in steps of two lattice sites, via the 
virtual occupation of the intermediate sites. 

Model - We consider a system of two electrons hop- 
ping on a chain of connected rhombi, as shown in FiglT^. 
When propagation from site to site is coherent this sys- 
tem can be modeled well by a single-band Hubbard 
model 



the Floquet theorem allows us to write solutions of the 



Hn = -J 



E 



J4'jkJ 



Cka + H.C. +C/^ "Jt'^j; . (1) 



where J is the intersite tunneling between nearest neigh- 
bors (j, k) , and U represents the energy cost of doubly- 
occupying a lattice site. The operators c]^ / c^^ are the 
usual creation / annihilation operators for an electron of 
spin cr on site j, and rijo- is the standard number op- 
erator. Each rhombus is threaded by a magnetic flux, 

giving rise to the AB phases exp [i4>jk\ on the tunnel- 
ing terms. In semiconductor quantum dots singlet-triplet 
mixing terms are typically rather weak, arising from the 
nuclear hyperfine interaction , and accordingly we ne- 
glect their effect here. We also confine our attention in 
this work to the singlet subspace, where one electron is 
spin-up and the other is spin-down, since in the triplet 
subspace the Hubbard-interaction is not operative. 

Results - In Fig[T}D we show the time-evolution of 
the system's wavefunction by plotting its overlap with 
the initial state, which consists of both electrons occu- 
pying the central site. In order to produce AB caging 
we set the applied flux to $ = tt In the absence of 
interactions {U = 0), the overlap displays regular sinu- 
soidal oscillations, indicating that the pair of electrons 
periodically revert to their original configuration. Exam- 
ining their dynamics in detail reveals that the electrons 
periodically oscillate from the initial site to its nearest 
neighbors, but propagate no further down the lattice due 
to the caging effect. Raising the interaction strength to 
U = O.IJ causes the AB-cage to partially open, allow- 
ing the electrons to spread over the whole lattice, and so 
causing the oscillations in the overlap to decay. As U is 
increased further this leakage occurs more rapidly. 

We now consider adding a periodically-oscillating po- 
tential that rises linearly along the lattice 



H{t) = Ho 



(2) 



Here E and uj parameterize the amplitude and frequency 
of the potential respectively, and Xj is the x-component of 
the location of site j. As H{t) is periodic in time, H{t) — 
H{t + T) where T = 27r/aj is the period of the driving, 



Schrodinger equation in the form 



Here e„ are the Floquet quasienergies, and \ipn{t)) are a 
set of T-periodic functions termed Floquet states. The 
quasienergies represent the appropriate generalization of 
energy eigenvalues to the case of a periodically-driven 
system, and have an analogous role in elucidating the sys- 
tem's behavior. In particular, exact or near-degeneracy 
of quasienergies results in the suppression of tunneling 
between the associated Floquet states: the phenomenon 
known as CDT [Ql]. 

High frequency (uj > U ) - In Figl2^ we show the 
Floquet quasienergies for a system in the high-frequency 
regime, where w is the dominant energy scale. We see 
that the quasienergies fall into two minibands; the higher 
miniband corresponding to Floquet states in which sites 
are doubly occupied, and the lower to states in which 
the electrons are separated in space. Both bands "col- 
lapse" at specific values of E/uj. In the high-frequency 
regime, a perturbative calculation of the Floquet system 
[lOl Il3| reveals that the driving potential has the effect 
of renormalizing the tunneling to an effective value J — > 
JqS — J Jo{Ex/u}), where x = 1/2 is the x-component 
of the spacing between neighboring lattice sites in units 
of the lattice spacing. Thus when E/uj 4.80, 11.04 .. ., 
which correspond to zeros of j7o, the effective tunneling 
is suppressed and the miniband collapses. 



The driving field thus provides a handle to directly 
control the ratio U/J^g, which can be enhanced by tun- 
ing the driving parameters near to a zero of JTq. This 
has been used, for example, to induce the Mott transi- 
tion in cold atom systems [IJ]. Another consequence of 
enhancing this ratio is the creation of repulsively-bound 
15[, or "doublons". Since binding is usually 
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associated with attractive forces, it may seem counter- 
intuitive for strong repulsive interactions to also produce 
this effect. Qualitatively it can be understood from an 
energetics argument. The energies of a one dimensional 
lattice form a Bloch band with a width of 2 J, and thus 
the maximum kinetic energy carried by two free parti- 
cles is 4 J. If the particles are initially prepared in a state 
with a potential energy much greater than 4 J, the initial 
state then cannot decay without the mediation of dis- 
sipative processes. In sufficiently clean systems the only 
decay path |ll| is via high-order scattering processes with 
unbound particles that carry off the excess energy, result- 
ing in a doublon lifetime that depends exponentially on 

U/JeS- 

To observe the dynamics of doublons we require Jeff 
to be close to zero (to enhance the ratio U/J^s for dou- 
blon formation to occur), but to be sufficiently large for 
the doublon tunneling to be no n- negligible. In Figl^b 
we show the occupation of the central site and its near- 
est neighbors for driving parameters E/uj — 4.80. For 
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FIG. 2: (a) Floquet quasienergies for a periodically-driven 
two electron system. Parameters are U = 2J, uj = 16J (the 
high-frequency regime). The quasienergies fall into two mini- 
bands: the lower (black) corresponds to Floquet states in 
which the electrons are separated, the higher (red) to states 
in which both electrons occupy the same site. The minibands 
are separated by the Hubbard gap U. Both minibands col- 
lapse at E/uj = 4.80, 11.04, corresponding to CDT. (b) For 
E/uj = 4.80 and $ = 7r/2 the occupation of the central site 
(solid black line) oscillates as electrons periodically tunnel to 
and from their nearest neighbors (see Fig[lji). For ease of 
comparison, the occupation of the neighboring site (dashed 
red line) is multiplied by 2. When $ is tuned away from 
this value AB-caging does not occur; the dot-dashed blue line 
shows the occupation of the central site for $ = for which 
the electrons diffuse throughout the lattice. 



general values of the magnetic flux $ the doublon slowly 
diffuses through the lattice, the two electrons remaining 
bound together, but not being localized. For the case of 
$ — 11/2^ however, the dynamics again shows a regular 
oscillatory motion due to AB-caging. Thus even though 
the system is strongly interacting [U — 2J) AB-caging 
can nonetheless be induced. Since the doublon has a 
charge of 2e, however, the caging now occurs when $ is 
equal to half a flux quantum. 

Resonant frequency - We now consider raising the in- 
teraction strength further. When w is no longer the dom- 
inant energy scale, the simple perturbation theory must 



El. This re- 



be generalized to include interactions 
veals that in this regime the effect of the driving field 
is particularly strong when it is resonant with the driv- 
ing frequency U = nuj, where n is an integer. Away 
from these specific frequencies the driving produces little 
structure in the quasienergy spectrum. 

In Fig[3^ we show the quasienergies obtained for pa- 
rameters U = 16J and uj = 16J: the n = 1 resonance. 
As before the Floquet states in which the electrons are 
separated fall into a miniband modulated by J'q{E/2uj). 
As predicted by perturbation theory |12l] . however, the 
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FIG. 3: (a) Resonant driving, U = 16J, lo = 16 J. Quasiener- 
gies corresponding to states in which the electrons are sep- 
arated (black) are modulated by the zeroth Bessel func- 
tion J'o{E /2ij). The quasienergies corresponding to doubly- 
occupied states (red) are instead modulated by J'i{E/2lj), 
and show band collapses at Ex/lj ~ 7.66, 15.44. (b) For 
E/uj = 7.50 and $ — 7r/2 AB-caging occurs, and the occu- 
pation of the central site (solid black line) oscillates from 2 
to zero, while correspondingly the occupation of the neigh- 
boring sites (dashed red line) oscillates between zero and 0.5. 
Tuning "1? away from this value destroys the AB-caging; the 
dot-dashed blue line shows the occupation of the central site 
for $ = 0. 



quasienergies in which the Hubbard interaction is opera- 
tive are instead modulated by Jn{E/2ui) where n is the 
order of the resonance. Consequently doublon formation 
now occurs when E I2uj is tuned near to a zero of jTi.. This 
is shown in Fig[3)D; when E/lu — 7.60, close to the first 
zero of Ji a stable doublon is formed, which undergoes 
AB-caging when the applied flux is set to = 7r/2. As 
before AB-caging only occurs for this speciflc flux value, 
and for other values of $ the doublon simply spreads 
through the lattice. It is interesting to note that un- 
like the high-frequency case, the oscillation is now com- 
plete; the central site empties completely as the dou- 
blon tunnels to the other sites in the AB-cage. This is 
a consequence of the resonance condition. In an effect 
analogous to photon assisted tunneUng, the doublon can 
absorb energy from the driving field to exactly compen- 
sate for [/, meaning that the system becomes effectively 
non-interacting. 



Low frequency (lo < U ) - Finally we consider the 
case when the driving frequency is much smaller than 
the interaction strength. If U/ J is sufficiently large then 
doublons will form in the static system, and the driving 
field can now be used to control their motion. In Figd^ 
we show the Floquet spectrum for states in which sites 
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FIG. 4: (a) For low driving frequencies, U — 16J, uj = 2 J, the 
quasienergies corresponding to doubly-occupied states form a 
miniband modulated by JoiE/uj). (b) Doublon dynamics for 
a 19-site system, where E/lj = 2.30 to suppress the direct 
doublon-tunneling. The doublon is initialized on the central 
site 1 (see below) and propagates along the spine of the lattice. 
After reaching site 3 (3') it is reflected, and returns to the 
center, (c) When the doublon tunneling is suppressed, motion 
occurs via second order tunneling; the doublon first unbinds 
(1) and then recombines (2). The doublon thus propagates in 
steps of two lattice sites, occupying the intermediate sites only 
virtually, (d) Schematic arrangement of the 19-site lattice. 
Filled circles indicate "spinal" sites, empty circles the "edge" 
sites. The second order tunneling clearly does not connect 
spinal sites to edge sites. 

are doubly-occupied, which we can regard as a miniband 
of doublon states. This miniband is clearly modulated by 
J'o{2Ex/uj), where the factor of 2 in the argument of the 
Bessel function occurs because the doublon has charge 2e. 
The miniband structure persists until E/uj > 8, which 
marks the onset of resonant driving [12]. For these pa- 
rameters this will be an n = 8 resonance (since U/lu — 8), 
but due to the strong driving potential required to reach 
this regime, we do not study it further here. 

Tuning E/uj close to 2.40 thus has the effect of sup- 
pressing direct tunneling of the doublon. Motion is still 
possible, however, via second order tunneling processes; 
the bound electrons hop to intermediate lattice sites (in- 
curring an energy cost of ?7), and then recombine after 
another hopping process, as illustrated schematically in 
FigHt. This motion has the unusual feature that it oc- 
curs in steps of two lattice sites, with the intermediate 
sites being occupied only virtually. If we separate the 
rhombus-chain lattice into "edge" and "spinal" sites, as 
shown in Fig|3Ji, it can be easily seen that this tunneling 
process does not connect spinal sites to edge sites, and 
so produces a separation between these two sub-lattices. 



In Figl^b we show the result of initializing the system 
in a doublon state at the center of a 19-site lattice. We 
can see that the doublon propagates smoothly along the 
spine of the lattice, the occupation of edge sites always 



remaining less than 0.006. This propagation occurs sym- 
metrically to the left and right, and so behaves as an 
electronic beam-splitter On reaching the final site 

(site 3, as shown in Fig|3ti) the wavepacket is reflected 
and revives in the central site with excellent fidelity. We 
note that the doublon does not reach site 4, although in 
principle the second-order tunneling process should per- 
mit this. This is essentially a finite-size effect ^3,]; the 
states that project onto the terminating sites (4 and 4') 
have a different symmetry to the other spinal states, and 
so lie in a different miniband, thereby isolating them from 
the dynamics of the spinal states. 

Conclusions - Combining AB-phases with a time de- 
pendent driving potential gives a extremely rich behav- 
ior. In particular, we have shown that AB-trapping can 
occur in an interacting system, by using CDT to convert 
pairs of electrons into doublons. The resonant behavior 
displayed by the driven system also provides a convenient 
means to measure the strength of interactions by suitably 
tuning the value of uj. Reducing the frequency of the 
driving gives complete coherent control over the dynam- 
ics of the doublons, allowing them to be localized within 
an AB-cage, or to propagate via an unusual second-order 
tunneling process. This permits the creation and control 
of spatially-separated entangled states of two electrons 
via the beam-splitter effect, with many potential appli- 
cations to quantum information. The fine discrimination 
of this phase-control method also allows doublons and un- 
paired electrons to be separately controlled, allowing in 
principle the coherent spatial separation of these compo- 
nents, and opening up exciting possibilities for quantum 
transport in interacting fcrmion systems. 
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